Phenomenological study of spin-triplet ferromagnetic 

superconductors 

Diana V. Shopova, Tsvetomir E. Tsvetkov, and Dimo I. Uzunov, 

CP Laboratory, Institute of Solid State Physics, 
<^> Bulgarian Academy of Sciences, BG-1784 Sofia, Bulgaria. 

O 

<N , 

^ . Abstract 

o 

Unconventional superconductivity with spin-triplet Cooper pairing is reviewed 
, on the basis of the quasi-phenomenological Ginzburg-Landau theory. The super- 

conductivity, in particular, the mixed phase of coexistence of ferromagnetism and 
unconventional superconductivity is triggered by the spontaneous magnetization. 
The mixed phase is stable whereas the other superconducting phases that usu- 



o 



ally exist in unconventional superconductors are either unstable, or, for particular 
values of the parameters of the theory, some of these phases are metastable at 
relatively low temperatures in a quite narrow domain of the phase diagram. The 
phase transitions from the normal phase to the phase of coexistence is of first 
order while the phase transition from the ferromagnetic phase to the coexistence 
phase can be either of first or second order depending on the concrete substance. 
\ The Cooper pair and crystal anisotropy are relevant to a more precise outline of 

the phase diagram shape and reduce the degeneration of the ground states of the 
system but they do not drastically influence the phase stability domains and the 
thermodynamic properties of the respective phases. 
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1 Inroduction 

The formation of Cooper pairs with a nonzero angular momentum was theoretically 
predicted [Tj in 1959 as a mechanism of superfluidity of Fermi liquids. In 1972 the 
same phenomenon - unconventional superfluidity due to a p-wave (spin triplet) Cooper 
pairing of 3 He atoms, was experimentally discovered in the mK range of temperatures; 
for details and theoretical description, see Refs. pOIHllIj- In contrast to the standard 
s-wave pairing in usual (conventional) superconductors, where the electron pairs are 
formed by an attractive electron-electron interaction due to a virtual phonon exchange, 
the widely accepted mechanism of the Cooper pairing in superfluid 3 He is based on 
an attractive interaction between the fermions ( 3 He atoms) due to a virtual exchange 
of spin fluctuations. Certain spin fluctuation mechanisms of unconventional Cooper 
pairing of electrons were proposed also for the depiction of discovered in 1979 heavy 
fermion superconductors (see, e.g., Refs. jSl El E|) as well as for some classes of high- 
temperature superconductors (see, e.g., Refs. [HI El Ell HH H21 HS1 CH HS1 HE! ) 
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The possible superconducting phases in unconventional superconductors are described 
in the framework of the general Ginzburg-Landau (GL) effective free energy func- 
tional PH] with the help of the symmetry groups theory. A variety of possible supercon- 
ducting orderings were predicted for different crystal structures [T7J HH1 HH I2D1 122 122] ■ 
A detailed thermodynamic analysis JUEIj °f the homogeneous (Meissner) phases and 
a renormalization group investigation ^1] of the superconducting phase transition up 
to the two-loop approximation have been also performed (for a three-loop renormal- 
ization group analysis, see Ref. [2S]; for effects of magnetic fluctuations and disorder, 
see PIES!)- 

In 2000, experiments [2E! at low temperatures (T ~ 1 K) and high pressure (T ~ 1 
GPa) demonstrated the existence of spin triplet superconducting states in the metallic 
compound UGe 2 . This superconductivity is triggered by the spontaneous magnetiza- 
tion of the ferromagnetic phase which exists at much higher temperatures and coexists 
with the superconducting phase in the whole domain of existence of the latter be- 
low T ~ 1 K; see also experiments published in Refs. [23 I2H1, and the discussion in 
Ref. - Moreover, the same phenomenon of existence of superconductivity at low 
temperatures and high pressure in the domain of the (T, P) phase diagram where the 
ferromagnetic order is present has been observed in other ferromagnetic metallic com- 
pounds (ZrZn 2 [HU] and URhGe pH] ) soon after the discovery [2H] of superconductivity 
in UGe 2 . 

In contrast to other superconducting materials, for example, ternary and Chevrel phase 
compounds, where the effects of magnetic order on superconductivity are also substan- 
tial (see, e.g., |H21 ESI EH E5j), in these ferromagnetic compounds the phase transition 
temperature (Tf) to the ferromagnetic state is much higher than the phase transition 
temperature (Tps) from ferromagnetic to a (mixed) state of coexistence of ferromag- 
netism and superconductivity. For example, in UGe2 Tp$ is 0.8 K whereas the critical 
temperature of the phase transition from paramagnetic to ferromagnetic state in the 
same material is Tf = 35K j2El I2Z| • It can be supposed that in such substances the 
material parameter T s defined as the (usual) critical temperature of the second order 
phase transition from normal to uniform (Meissner) supercondicting state in a zero ex- 
ternal magnetic field is much lower than the phase transition temperature Tps- Note, 
that the mentioned experiments with the compounds UGe 2 , URhGe, and ZrZn 2 do 
not give any evidence for the existence of a standard normal-to-superconducting phase 
transition in a zero external magnetic field. 

Moreover, it seems that the superconductivity in the metallic compounds, mentioned 
above, always coexists with the ferromagnetic order and is enhanced by it. As claimed 
in Ref. |2E] in these systems the superconductivity seems to arise from the same elec- 
trons that create the band magnetism, and is most naturally understood as a triplet 
rather than spin-singlet pairing phenomenon. Note, that all three metallic compounds, 
mentioned so far, are itinerant ferromagnets. A similar type of unconventional su- 
perconductivity has been suggested [US] as a possible outcome of recent experiments 
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in Fe [37], where a superconducting phase is discovered at temperatures below 2 K 
for pressures between 15 and 30 GPa. Note, that both vortex and Meissner super- 
conductivity phases [SZj have been found in the high-pressure crystal modification of 
Fe which has a hexagonal close-packed crystal structure. In this hexagonal lattice the 
strong ferromagnetism of the usual bcc iron crystal probably disappears jSHl- Thus one 
can hardly claim that there is a coexistence of ferromagnetism and superconductivity 
in Fe but the clear evidence of a superconductivity is also a remarkable achievement. 

The important point in all discussions of the interplay between superconductivity and 
ferromagnetism is that a small amount of magnetic impurities can destroy superconduc- 
tivity in conventional (s-wave) superconductors by breaking up the (s-wave) electron 
pairs with opposite spins (paramagnetic impurity effect [EE!)- I* 1 this aspect the phe- 
nomenological arguments [EH] and the conclusions on the basis of the microscopic theory 
of magnetic impurities in s-wave superconductors [38J are in a complete agreement with 
each other; see, e.g., Refs. jE2 ESI EH EH] • In fact, a total suppression of conventional 
(s-wave) superconductivity should occur in the presence of an uniform spontaneous 
magnetization M, i.e., in a standard ferromagnetic phase [HE]- The physical reason for 
this suppression is the same as in the case of magnetic impurities, namely, the opposite 
electron spins in the s-wave Cooper pair turn along the vector M in order to lower 
their Zeeman energy and, hence, the pairs break down. Therefore, the ferromagnetic 
order can hardly coexist with conventional superconducting states. Especially this is 
so for the coexistence of uniform superconducting and ferromagnetic states when the 
superconducting order parameter ip(x) and the magnetization M(x) do not depend 
on the spatial vector x. 

But yet a coexistence of s-wave superconductivity and ferromagnetism may appear 
in uncommon materials and under quite special circumstances. Furthermore, let us 
emphasize that the conditions for the coexistence of nonuniform ("vertex", "spiral", 
"spin-sinosoidal" or "helical") superconducting and ferromagnetic states are less re- 
strictive than those for the coexistence of uniform superconducting and ferromagnetic 
orders. Coexistence of nonuniform phases has been discussed in details, in partic- 
ular, experiment and theory of ternary and Chevrel-phase compounds, where such 
coexistence seems quite likely; for a comprehensive review, see, for example, Refs. 

[321 E31 Ell ESI E0]. 

In fact, the only two superconducting systems for which the experimental data allow 
assumptions in favor of a coexistence of superconductivity and ferromagnetism are 
the rare earth ternary boride compound ErRh 4 B 4 and the Chervel phase compound 
HoMo 6 S 8 ; for a more extended review, see Refs. [EE1EI]- In these compounds the phase 
of coexistence most likely appears in a very narrow temperature region just below the 
Curie temperature Tf of the ferromagnetic phase transition. At lower temperatures 
the magnetic moments of the rare earth 4/ electrons become better aligned, the mag- 
netization increases and the s-wave superconductivity pairs formed by the conduction 
electrons disintegrate. 
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We shall not extend our consideration over all important aspects of the long stand- 
ing problem of coexistence of superconductivity and ferromagnetism rather we shall 
concentrate our attention on the description of the newly discovered coexistence of 
ferromagnetism and unconventional (spin-triplet) superconductivity in the itinerant 
ferromagnets UGe 2 , ZrZn 2 , and URhGe. Here we wish to emphasize that the main 
object of our discussion is the superconductivity of these compounds and, in a second 
place in rate of importance we put the problem of coexistence. The reason is that 
the existence of superconductivity in such itinerant ferromagnets is a highly nontrivial 
phenomenon. As noted in Ref. [12] the superconductivity in these materials seems 
difficult to explain in terms of previous theories [H21 133 ESj and requires new concepts 
for the interpretation of experimental data. 

We have already mentioned that in ternary compounds the ferromagnetism comes from 
the localized 4/ electrons whereas the s-wave Cooper pairs are formed by conduction 
electrons. In UGe2 and URhGe the 5/ electrons of U atoms form both superconducting 
and ferromagnetic order [2H1 EH] • In ZrZn 2 the same twofold role is played by the 4d 
electrons of Zr. Therefore the task is to describe this behavior of the band electrons at 
a microscopic level. One may speculate about a spin-fluctuation mediated unconven- 
tional Cooper pairing as is in case of 3 He and heavy fermion superconductors. These 
important issues have not yet a reliable answer and for this reason we shall confine our 
consideration to a phenomenological level. 

Reliable experimental data, for example, the data about the coherence length and 
the superconducting gap |2H1 I2H Ell EH] , are in favor of the conclusion about a spin- 
triplet Cooper pairing in these metallic compounds, although the mechanism of the 
pairing remains unclear. We shall essentially use this reliable conclusion. Besides, 
this point of view is consistent with the experimental observation of coexistence of 
superconductivity only in the low temperature part of the ferromagnetic domain of the 
phase diagram (T, P)which means that a pure (non ferromagnetic) superconducting 
phase has not been observed. This circumstance is also in favor of the assumption of 
a spin-triplet superconductivity. Our investigation leads to results which confirm this 
general picture. 

On the basis of the experimental data and conclusions presented for the first time 
in Refs. |2T)] and shortly afterwards confirmed in Refs. (23 I2H1 EDI EH] it can be 
accepted that the superconductivity in these magnetic compounds is considerably en- 
hanced by the ferromagnetic order parameter M and, perhaps, it could not exist 
without this "mechanism of ferromagnetic trigger," or, in short, "TVf -trigger." Such 
phenomenon is possible for spin-triplet Cooper pairs, where the electron spins point 
parallel to each other and their turn along the vector of the spontaneous magnetization 
M does not cause a break down of the spin-triplet Cooper pairs but rather stabilizes 
them and, perhaps, stimulates their creation. We shall describe this phenomenon at a 
phenomenological level. 

Recently, the phenomenological theory that explains the coexistence of ferromagnetism 
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and unconventional spin-triplet superconductivity of GL type was developed [121 E3] • 
The possible low-order couplings between the superconducting and ferromagnetic order 
parameters were derived with the help of general symmetry group arguments and 
several important features of the superconducting vortex state in the ferromagnetic 
phase of unconventional ferromagnetic superconductors were established |%2"1 . 

In this paper we shall use the approach presented in Refs. (13 EH] to investigate the 
conditions for the occurrence of the Meissner phase and to demonstrate that the pres- 
ence of ferromagnetic order enhances the p-wave superconductivity. Our consideration 
is focused on the ground state, namely, we are interested in uniform phases, where 
the order parameters (the superconducting order parameter ip and the magnetization 
vector M = {Mj,j = 1,2,3}), do not depend on the spatial vector x e V, where V 
is the volume of the system. Recent results about the phase diagram and the phase 
transitions jUJ ES] , and thermodynamic quantities [IE] will be essentially used in our 
investigation. 

Our study is based on the mean-field approximation ^3] as well as on familiar results 
for the possible phases in nonmagnetic superconductors with triplet (p-wave) Cooper 
pairs [HI H21 HH|- Results from Refs.jSl EE] EE] Ej will be reviewed and extended. 
In our preceding investigation |HJ EE1 EE] both Cooper pair anisotropy and crystal 
anisotropy have been neglected in order to clarify the main effect of the coupling be- 
tween the ferromagnetic and superconducting order parameters. The phenomenological 
GL free energy is quite complex and the inclusion of these anisotropics is related with 
lengthy formulae and a multivariant analysis which obscures the final results. Here we 
shall point our attention to the effect of the Cooper pairs anisotropy. 

There exists a formal similarity between the phase diagram we have obtained and the 
phase diagram of certain improper ferroelectrics [JH1 EE EDI OH] E2 IHSj- The variants 
of the theory of improper ferroelectrics, known before 1980, were criticized in Ref. jSHl 
for their oversimplification and inconsistency with the experimental results. But the 
further development of the theory has no such disadvantages (see, e.g., Ref. [SUES)- 
We should emphasize that the symmetry of the GL model of spin-triplet ferromagnetic 
superconductors is quite different from the symmetry of known models in ferroelectrics 
and, hence, the results for ferroelectric systems can hardly be applied to superconduc- 
tors without additional investigations. 

2 Ginzburg-Landau free energy 

Consider the GL free energy F(ip,M) = Vf(ip,M) , where the free energy density 
f(ip,M) (for short hereafter called "free energy") of a spin-triplet ferromagnetic su- 
perconductor is given by 
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= a s H 2 + ||^| 4 + y|^| 2 + |E^| 4 + a / M2 + y M4 C 1 ) 

j'=i 

In Eq. (1), tj) — {ipj]j = 1, 2, 3} is a three-dimensional complex vector = ip'- + iip") 
describing the unconventional (spin-triplet) superconducting order and B = (H + 
AttM) = V x A is the magnetic induction; H = {Hj;j = 1,2,3} is the external 
magnetic field, A = {A,; j = 1,2,3} is the magnetic vector potential (V. A = 0). In 
Eq. (1), b s > 0, bf > 0, a/ = aif(T — Tf) is given by the positive material parameter af 
and the ferromagnetic critical temperature Tj corresponding to a simple ferromagnet 
(ip = 0), and a s = a s (T — T s ), where a s is another positive material parameter and 
T s is the critical temperature of a standard second order phase transition which may 
occur at \H\ = M. = 0; M. = \M\. The parameter u s describes the anisotropy of the 
spin-triplet Cooper pair whereas the crystal anisotropy is described by the parameter 

v s dUEEEl. 

The two orders - the magnetization vector M = {Mj} and ip = {Aj}, interact through 
the last two terms in (1). The 70— term ensures the triggering of the supercon- 
ductivity by the ferromagnetic order (70 > 0) whereas the 5— term makes the model 
more realistic in the strong coupling limit 02]. Both ^M-interaction terms included 
in (1) are important for a correct description of the temperature-pressure (T, P) phase 
diagram of the ferromagnetic superconductor [¥3J 03] . In general, the parameter 5 for 
ferromagnetic superconductors may take both positive and negative values. 

As we are interested in the ground state properties, we set the external magnetic 
field equal to zero (H = 0). Besides, we emphasize that the magnetization vector M 
may produce vortex superconducting phase in case of type II superconductivity. The 
investigation of nonuniform (vortex) states can be made with the help of gradient terms 
in the free energy which take into account the spatial variations of the order parameter 
field ip. This task is beyond our present consideration. Rather we investigate the basic 
problem about the possible stable uniform (Meissner) superconducting phases which 
may coexist with uniform ferromagnetic order. For this aim the free energy (1) is quite 
convenient. 

In case of a strong easy axis type of magnetic anisotropy, as is in UGe2 j2E], the 
overall complexity of mean-field analysis of the free energy f(ip,M) can be avoided 
by performing an "Ising-like" description: M = (0,0,.M). Further, because of the 
equivalence of the "up" and "down" physical states (±M) the thermodynamic analysis 
can be done within the "gauge" M. > 0. But this stage of consideration can also be 
achieved without the help of crystal anisotropy arguments. When the magnetic order 
has a continuous symmetry one may take advantage of the symmetry of the total free 
energy f(ip,M) and avoid the study of equivalent thermodynamic states that occur 
as a result of the respective symmetry breaking at the phase transition point but 
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they have no effect on thermodynamics of the system. In the isotropic system one may 
again choose a gauge, in which the magnetization vector has the same direction as z-axis 
(\M\ = M z = M) and this will not influence the generality of thermodynamic analysis. 
Here we shall prefer the alternative description within which the ferromagnetic state 
may occur as two equivalent "up" and "down" domains with magnetizations M. and 
—A4, respectively. 

We shall use adequate notations which reduce the number of parameters. With the 
help of 

b = (b s + u s + v a ) > (2) 
we redefine the order parameters and the other parameters in the following way: 

!fj = = fae 9 * , M = b) /A M , (3) 

a s a f u s v s 

r= V5' t= ^V,' U ' = 1T' V = T' 

- 7 ° S 

7 ~ bv%f ' 71 ~ (bb f y/ 2 ■ 

Having in mind our approximation of uniform ip and M and the notations (2) - (3), 
the free energy density f(ip, M) can be written in the form 

f{ip,M) = r0 2 + ^0 4 + 270 1 2 Msin(^-^ 1 )+ 7l 2 M 2 + tM 2 + iM 4 (4) 
-2w [0 2 2 sin 2 (0 2 - 9 X ) + <j>l<f%an 2 (9x - 6 3 ) + 2 2 sin 2 (^ 2 - 9 3 )] 

In the above free energy the order parameters ip and M are defined per unit volume. 

We assume that Tf > T s . This is the case when the superconductivity is triggered 
by the magnetic order. We shall discuss the stable phases in the temperature region 
T > T s . The case Tf < T s may be also important for neutron stars so it needs a special 
investigation. When T s ~ Tf a quite simple analytical treatment is possible. All these 
cases may be of interest to the description of ferromagnetic superconductivity in stellar 
objects whereas in condensed matter only case of Tf ^> T s has been observed so far. 

We work in the framework of the standard mean-field analysis jT^j. The stable phases 
correspond to global minima of the GL energy (1). The equilibrium phase transition 
line separating two phases is defined by the thermodynamic states, where the respective 
GL free energies are equal. 

3 Phases 

We shall not enumerate and discuss all phases described by Eq. (1). Rather we shall 
focus our attention on the stable phases at relatively high temperatures (T > T s ). The 
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calculations show that for temperatures T > T s , i.e., for r > 0, there are three stable 
phases. Two of them are quite simple: the normal phase (ip = M = 0) with existence 
and stability domains given by t > and r > 0, and the ferromagnetic phase (FM) 
given by ip = and M 2 = —t whose existence condition is t < 0,the is stability domain 
defined by the inequalities r > jit and 



r > 7jt + 7V-t • (5) 

The third stable phase is the phase of coexistence of superconductivity and ferromag- 
netism (hereafter referred to as FS). This phase is the main object of our consideration. 
It is given by the following equations: 

01 = 4>2 = ~J= , 03 = , (6) 

<f) 2 = (± 7 M - r - 7l M 2 ) , (7) 
(1 - 7l 2 )M 3 ± ^77iM 2 +(t--L- ^r) M ± ^ = , (8) 



and 

(# 2 -#i) = ^ + 27rA;, (9) 

(fc = 0,±1,...). The upper sign in Eqs. (7) - (9) corresponds to a domain where 
sin (#2 — 9\) = —1 and the lower sign corresponds to a second domain which may be 
referred to as FS*; there sin(# 2 — 0\) = 1. These two domains are equivalent and 
describe the same ordering. We shall focus on the upper sign in (7) - (9), i.e. on FS. 
The phase diagram (t, r) has qualitatively the same shape as the phase transition lines 
corresponding to w = but there are essential quantitative differences between them. 
We shall discuss them in the next section. Note, that the system exhibits both first 
and second order phase transitions and complex phase transition points: triple and 
tricritical points |Hl HZ] . 



4 Anisotropy effects 

Our analysis demonstrates that when the anisotropy of the Cooper pairs is taken into 
account, there will be not drastic changes in the shape the phase diagram for r > 
and the order of the respective phase transitions. Of course, there will be some changes 
in the size of the phase domains and the formulae for the thermodynamic quantities. 
Besides, and this seems to be the main anisotropy effect, the w- and v -terms in the free 
energy lead to a stabilization of the order along the main crystal directions which, in 
other words, means that the degeneration of the possible ground states is considerably 
reduced. So there will be a smaller number of marginally stable states. 
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Let us neglect the crystal anisotropy by setting v s = in Eq. (1) and concentrate our 
attention on the Cooper pair anisotropy described by the w s -term in the GL model. 
The dimensionless anisotropy parameter w ~ u s given by Eq. (3) can be either positive 
or negative depending on the sign of u s . Obviously when u s > 0, the parameter w will 
be positive too (0 < w < 1). We shall illustrate the influence of Cooper-pair anisotropy 
in this case. The order parameters (M, <pj, 6j) are given by Eqs. (6), (9), 



± 1 M - r - 7iM 2 



w 



and 



[l-w- 7l 2 )M 3 ± -771M 2 + 



> 



r 



t(l -w)- — -^ x r 







(10) 



11^ 



where the meaning of the upper and lower sign is the same as explained just below 
Eq. (9). We consider the FS domain corresponding to the upper sign in the Eq. (10) 
and (11). The stability conditions for FS read, 



(2 - w)-fM - r - 7xM 2 
1 — w 



> 



(12) 



and 



W 



7M -wr- wjiM 2 > 



3(1 - w - 7i)M 2 + 3771M + t(l -w)- 



r 



7ir 



> 0. 



(13) 



(14) 



For M 7^ (7/271) we can express the function r(M,t) defined by Eq. (11), substitute 
the obtained expression for r(M,t) in the existence (10) and stability conditions (12)- 
(14) and do the analysis in the same way as for w = 47 . The most substantial 
qualitative difference between the cases w > and w < is that for w < the stability 
of FS is bounded for r < 0. 



5 Conclusion 

We have done an investigation of the M-trigger effect in unconventional ferromagnetic 
superconductors. This effect due to the Mipi ^-coupling term in the GL free energy 
consists of bringing into existence of superconductivity in a domain of the phase dia- 
gram of the system that is entirely in the region of existence of the ferromagnetic phase. 
This form of coexistence of unconventional superconductivity and ferromagnetic order 
is possible for temperatures above and below the critical temperature T s , which corre- 
sponds to the standard phase transition of second order from normal to Meissner phase 
- usual uniform superconductivity in a zero external magnetic field, which appears out- 
side the domain of existence of the ferromagnetic order. Our investigation has been 
mainly intended to clarify the thermodynamic behaviour at temperatures T s < T < Tf, 
where the superconductivity cannot appear without the mechanism of M-triggering. 
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We have described the possible ordered phases (FM and FS) in this most interesting 
temperature interval. 

The Cooper pair and crystal anisotropies have also been investigated and their main 
effects on the thermodynamics of the triggered phase of coexistence are established. 
In discussions of concrete real materials one should take into account the respective 
crystal symmetry but the variation of the essential thermodynamic properties with the 
change of the type of symmetry is not substantial when the low symmetry and low 
order (in both M and if)) 7-term is present in the free energy. 

Below the superconducting critical temperature T s a variety of pure superconducting 
and mixed phases of coexistence of superconductivity and ferromagnetism exists and 
the thermodynamic behavior at these relatively low temperatures is more complex 
than in known cases of improper ferroelectrics. The case Tf < T s also needs a special 
investigation. Our results are referred to the possible uniform superconducting and 
ferromagnetic states. Vortex and other nonuniform phases need a separate study. 

The relation of the present investigation to properties of real ferromagnetic compounds, 
such as UGe2, URhGe, and ZrZn 2 , has been discussed throughout the text. In these 
real compounds the ferromagnetic critical temperature is much larger than the super- 
conducting critical temperature (Tf ^> T s ) and that is why the M-triggering of the 
spin-triplet superconductivity is very strong. Moreover, the 71-term is important to 
stabilize the FM order up to the absolute zero (0 K), as is in the known spin-triplet fer- 
romagnetic superconductors. The neglecting [13] of the symmetry conserving 7 1 -term 
hinders the proper description of real substances of this type. More experimental infor- 
mation about the values of the material parameters (a s ,df, ...) is required in order to 
outline the thermodynamic behavior and the phase diagram in terms of thermodynamic 
parameters T and P. In particular, a reliable knowledge about the dependence of the 
parameters a s and aj on the pressure P, the value of the characteristic temperature T s 
and the ratio a s /aj at zero temperature are of primary interest. 
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